Theoretical search for superconductivity in SC3XB perovskites and weak 
ferromagnetism in Sc^X (X = Tl, In, Ga, Al). 
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A possibility for a new family of intermetallic perovskite superconductors ScsXB, with 
X = Tl, In, Ga and Al, is presented as a result of KKR electronic structure and pseudopoten- 
tial phonon calculations. The large values of computed McMillan-Hopfield parameters on scandium 
suggest appearance of superconductivity in Sc^XB compounds. On the other hand, the possibil- 
ity of weak itinerant ferromagnetic behavior of ScaX systems is indicated by the small magnetic 
moment on Sc atoms in two cases of X = Tl and In. Also the electronic structure and resulting 
superconducting parameters for more realistic case of boron-deficient systems ScsXBx are com- 
puted using KKR-CPA method, by replacing boron atom with a vacancy. The comparison of the 
calculated McMillan-Hopfield parameters of the ScaXB series with corresponding values in MgCNis 
and YRhsB superconductors is given, finding the favorable trends for superconductivity. 



I. INTRODUCTION 

The motivation to search for a superconductivity in the 
intermetallic series of ScaXB compounds was prompted 
by the very interestiiig and quite non-typical supercon- 
ductivity of MgCNia [l| perovskite, with the critical tem- 
perature Tc ~ 8 K. First of all, the large amount of 
nickel atoms would suggest ferromagnetic properties in 
this compound and propensity to form magnetic state 
is really large, due to the van Hove singularity near the 
Fermi level This may be a reason of some unex- 

pected effects, which were observed experimentally in 
this material (t^I, despite MgCNja belongs 

to electron-phonon type superconductors 0, with a 
very large isotope effect on carbon . 

Among other intermetallic perovskites, similar to 
MgCNia, a single compound was earlier reported to be 
a superconductor: YRhsB with = 0.76 K How- 
ever, since discovery of superconductivity in MgCNis, no 
other superconductor has been found so far. Systematic 
experimental study showed some indications of su- 
perconductivity in CaBa;Pd3 (Tc ~ 1 K) and NbBa;Rh3 
(Tc ~ 6 K), but the superconducting phase was finally 
not identified. Also, no evidence of superconductivity 
was found experimentally in ZnCNia i despite the cal- 
culated electronic structure beeing very close to that of 
MgCNis. Here the carbon deficiency was proposed as a 
possible explanation [l3j| . 

The predictions of superconductivity in one of the 
compounds from the entitled series - ScalnB - was 
announced in our conference paper [l4j . Interesting 
properties of ScalnB encouraged us to investigate elec- 
tronic structure and electron-phonon coupling along the 
whole isoelectronic series of ScaXB compounds, with 
X = Tl, In, Ga, Al, and detailed results of our theo- 



retical study is presented in this work. 

The crystal structure of ScsATB (space group Pm-Srn, 
CaTiOs type) is given in Fig. n The synthesis of Tl- 
and In-containing cubic perovskites was reported |l5j . 
while Ga- and Al-containing compounds are still hypo- 
thetical. The perovskite structure of SC3XB may also 
be viewed as a CusAu-type cubic SC3X, with additional 
boron atom placed in the center (Fig. QJ. In particular, 
the Scsin system was found to exist, in two polymorphic 
forms 16]: one is the above mentioned, poorly known 
cubic CusAu, detected under a high pressure 117|, and 
the second one is hexagonal Ni^Sn-type. in which weak 
ferromagnetic properties were observed |l8l |. The ferro- 
magnetism of hexagonal Scain guided us to study pos- 
sibility of magnetic behaviors in the binary series of cu- 
bic ScaAT compounds. Indeed, the computations shows 
that the cubic Scain is magnetic as well as the hexagonal 
phase. For other systems the results of electronic struc- 
ture calculations are also presented, although we are not 
aware of the existence of these compounds. Noteworthy, 
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FIG. 1: Crystal structure of the ScsXB perovskite. The 
atomic positions are: X: la (0,0,0), B: lb (1/2,1/2,1/2) and 
Sc: 3c (0,1/2,1/2), (1/2,0,1/2), (1/2,1/2,0). If central B atom 
is removed, perovskite structure changes into CusAu-type. 
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spin-polarized computations of SC3XB shows the ground 
state to be non-magnetic. We also intend to underline 
some correlations between tendency to superconductiv- 
ity in SC3XB on one hand and to weak ferromagnetism 
in the corresponding Sc^X on the other hand. 



II. CALCULATIONS 

First-principle calculations of the superconducting pa- 
rameters, due to the subtle nature of such phenomenon, 
are still challenging problem in the Density-Functional 
computations. McMillan 19] in his fundamental work 
showed, that electron-phonon coupling (EPC) calcula- 
tions may be decoupled into electronic and phonon con- 
tributions, when using a few approximations. The com- 
monly used formula for the essential parameter for su- 
perconductivity - EPC constant A, resulting from gen- 
eralization of his concept to multi-atomic compounds, is 
as follows: 



A 



E 



(1) 



where i corresponds to all atoms in the unit cell, with rji 
being the electronic, and Mi{ujf) lattice contribution to 
the electron-phonon interaction parameter. 

The electronic part of the EPC constant - McMillan- 
Hopfield parameters rji 0, |23, which describe the re- 
sponse of electrons from Fermi surface to displacements 
of atoms, are computed within the Gaspari-Gyorffy "^ll 
method, using so-called Rigid-Muffin-Tin Approxima- 
tion (RMTA). The relevant formula for 

V is _22] (i sub- 
script is dropped): 



rj = 



E 



{2l + 2)ni{EF)ni+i{EF) 
{21 + 1){21 + 3)n{EF) 



ar 



(2) 

where V{r) is the self-consistent spherically-symmetric 
potential on given atom site, Rmt is the radius of the 
i-th muffin-tin sphere, Ri{r) is a regular solution of the 
radial Schrodinger equation (normalized to unity inside 
the MT sphere), and ni{EF) is the ?-th partial density 
of states at the Fermi level {Ep) on the site considered. 
This method usually involves three main assumptions 

(i) rigid-ion approximation, in which potential inside the 
MT sphere moves rigidly with the ion, and the change 
in crystal potential, caused by the atom displacement, is 
given by the potential gradient; 

(ii) local-vibration approximation - more generally A = 
^ij^ where {i,j) refer to two atoms in the unit cell; 

in this approximation off-diagonal terms are neglected. 



A,,; 



5ij A, 



(Eq. O) 



(iii) spherical hand approximation, which leads to only 
dipole transitions ^ / + 1) in Eq. 

It is known that all mentioned assumptions give sat- 
isfactory results for the transition metal elements and 



cubic-site symmetry. In simple metals the assumptions 
(i) and (ii) generally underestimate ||2J] EPC due to a 
poorly screened crystal potential, whereas, the correc- 
tions to (iii) are expected to be small for cubic transition 
metals |25l |. 

In order to allow for independent calculations of elec- 
tron and phonon subsystems, another simplification is 
necessary. Estimation of {iof) in Eq. 0J, using ex- 
clusively phonon density of states F{uo) instead of full 
electron-phonon coupling function a^F{uj), requires an 
assumption, that the electron-phonon interaction is con- 
sidered to be independent of phonon frequency u, so 
the electron-phonon interaction factor a^iuj) cancels [7ll | 
when calculating {ujf) (see e.g. Ref. H^). 

The assumption that EPC magnitude does not change 
with phonon frequency is well fulfilled e.g. in niobium 
|27j , but may not be of the same behavior in multi-atom 
compounds. However, the RMTA method was success- 
fully used for analyzing EPC in many superconducting 
materials, li ke p ure metals [l^, binary alloys [29| . ^-15 
compounds f3(l| , transition metal carbid es [3lL l3^ , boro- 
carbides jS^] or metal-hydrogen system [sj . Reasonable 
results (as far as A is concerned) were usually obtained, 
even in such unusual superconductor, as MgB2 This 
formalism was also helpful in discussing phonon-based ef- 
fects in high-temperature superconductors [2^ [s^ . Cer- 
tainly, in order to better understand the electron-phonon 
interactions in a superconductor, one has to perform 
more elaborated calculations, using the Eliashber g ga p 
equations (see, for example, results for MgB2 [stI l3Sj|). 
Nevertheless, such electronic structure and phonon calcu- 
lations, together with the simplified RMTA framework, 
are very useful and efficient tools of early stage of looking 
for new superconducting systems. 

In line with this methodology, McMillan-Hopfield pa- 
rameters and phonon frequency moments were calculated 
for ScaATB compounds, and then, using Eq. EPC pa- 
rameters A were deduced. 

Electronic structure calculations were performed with 
the KKR method, which in case of disordered systems 
(Sec. IIIB|I was implemented together with coherent po- 
tential approximation (CPA), as described in details in 
Ref. js^. Additionally, for ordered compounds com- 
putations were performed with the full-potential (FP) 
code, based on the FP-KKR formalism, widely discussed 
by many authors l4(]|. with technical details shown e.g. 
when applying to Si |4j|, or to the electric field gradient 
calculations |42| |. In our practice extracting bands was 
done with the novel quasi-linear algorithm which 
allows for more precise and less time-consuming band 
structure calculations, comparing to conventional tech- 
niques. The crystal potential was constructed within the 
local density approximation (LDA) and Perdew-Wang 
formula was used for the exchange-correlation part. 
For all calculations angular momentum cut-off Imax = 3 
was set. Highly converged results were obtained for ^120 
k-points grid in the irreducible part of Brillouin zone 
(IRBZ), but they were also checked for convergence us- 
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ing more dense k-mesh and Imax = 4. Electronic den- 
sities of states were computed using a tetrahedron k~ 
space integration technique (up to 700 tetrahedrons in 
IRBZ). In the case of investigated compounds, the dif- 
ferences between the full potential and spherical poten- 
tial KKR calculations of the density of states (DOS) 
values at the Fermi level were of order of a few per- 
cents. The justification for using different kinds of ap- 
proximations, when computing McMillan-Hopfield pa- 
rameters, was supported by checking the numerical val- 
ues on changing different input data (like nonrelativis- 
tic vs. semirclativistic calculations, type of exchange- 
correlation potential, radii of muffin-tin spheres), which 
we also found to differ on the range of a few percentage. 
Such accuracy we consider satisfactory to explain trends 
in the number of systems, but for more detailed analy- 
sis the fine computation with the full potential will be 
necessary, especially if the Fermi surface specific features 
become important. 

Phonon calculations were undertaken for realistic ap- 
proximation of the phonon part of the electron-phonon 
couphng constant A, i.e. juj'j) parameter in Eq. We 
used the P Wscf package [i^ , where the plane wa ve p seu- 
dopotential technique and perturbation theory |46l | was 
implemented. For Sc, B, and Tl atoms ultrasoft pseu- 
dopotentials were employed, for Al, Ga and In norm- 
conserving pseudopotentials were taken. LDA parame- 
terization of Perdew and Zunger |47j| was implemented. 
Plane-wave kinetic energy and charge density cut-offs 
were set to 30 Ry and 350 Ry, respectively. The Brillouin 
zone integration smearing technique of Methfessel and 
Paxton l48| (with parameter a = 0.02 Ry) was used dur- 
ing the calculations. In order to obtain the phonon DOS 
F{uj), first the dynamical matrices on (5,5,5) q-point grid 
were calculated. Then, using the Fourier transformation 
on the same grid, real-space interatomic force constants 
were computed. The final result - phonon densities of 
states - were obtained from frequencies calculated from 
force constants on (10,10,10) grid and using the tetrahe- 
dron method. The combined results of the KKR elec- 
tronic structure study and the phonon DOS calculations 
were then used to estimate the EPC constant A. 

We also estimate the superconducting transition tem- 
perature Tc from McMillan formula Eq. (O, using 
the modified factor (uj)/1.2 



TABLE I: Lattice parameters in ScsXB series, units; 
1 Bohr = 0.5292 A. 



M 

1.2 



exp 



1.04(1 + A) 



A-^*(1 + 0.62A) 



(3) 



The absolute value of Tc depends also on Coulomb 
pseudopotential parameter jj*, and its influence on Tc 
magnitude is additionally discussed. 



A. Predictions of superconductivity in SC3XB 

Crystal structure and atomic positions of the SC3AB 
system are presented in Fig. U For all investigated or- 



Compound 


a experimental fl5] 


ao calculated 


SC3TIB 


8.541 


8.650 


Sc3lnB 


8.618 


8.610 


ScaGaB 




8.550 


SC3AIB 




8.610 



dered SC3AB compounds, lattice constants were calcu- 
lated from total energy minimum, as a first step of the 
PWscf phonon calculations. Both KKR electronic struc- 
ture and phonon computations were then performed with 
these equilibrium lattice parameters. Except for X — Al, 
semi-relativistic calculations results are presented here. 
In order to obtain the McMillan-Hopfield parameters, 
the following muffin-tin spheres radii (in the lattice con- 
stant ao unit) were employed: Rsc = 0.325, Rx = 0.36 
and Rb = 0.174. Eq. (0) suggests, that ?]i parameter 
may be sensitive to a choice of Rmt, so we also checked 
the influence of the computational geometry on final re- 
sults. The Rmt spheres variation in the range of 15% 
changed the results in about 5%. The difference in ob- 
tained results for somehow arbitrary MT geometry will 
be treated as an computational accuracy of the rji value. 

In Table results of lattice constant optimization for 
SC3AB are summarized. In both existing compounds the 
agreement between experimental and theoretical values 
is quite good. The variation of ao with atom X reflects 
the change in ionic radius of the X element, where Ga 
has the smallest radius, and Tl the largest. For SC3TIB 
calculated value ao is larger than the experimental one, 
which is not so common case. Similar effect was also 
observed in YRh3B 51, 52] . The smaller experimental 
value of a in SC3TIB could indicate that the measured 
sample was boron deficient, as it was earlier suggested 
for YRhgB 52]. 

Figure 121 presents electronic DOS in the ScaXB com- 
pounds. As we can see, the DOS shape is very similar 
along this series, due to the same number of valence elec- 
trons. Noteworthy, even SC3AIB, where Al contains only 
s— an p— like orbitals, exhibits electronic structure quite 
similar to other SC3AB compounds, containing complete 
nd^'^ shell from X atom. This shell forms semicore level 
located about 1 Ry below Ep (not shown). It supports 
widely accepted conclusion, that 3d-Ga, 4c?-In and 5d-T\ 
electrons build core-like levels, and the effect on upper 
lying electronic states (particularly near Ep) is small. 

The most important part of the density of states near 
the Fermi level in all compounds is formed from Sc d- 
states hybridized with p-states from B and X. The 
site-contributions to n{Ep) are presented in Tab.HTI 7^. 
Electronic dispersion curves for the representative com- 
pound - ScsInB - are shown in Fig.|3| and they are very 
similar in other cases (not given). The DOS sharp peak 
below Ep comes from a flat band, which is best seen in 
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TABLE 11: Electronic and phonon properties of ScaXB. 
nlEp) is given in (1/Ry/atom), 77 in (mRy/Bohr^/atom), lo 
in (meV), M{cj^> in (mRy/Bohr^). 



-0.45 -0.3 -0.15 
Energy (Ry) 



FIG. 2: KKR total DOS of isoelectronic ScsXB perovskites 
calculated at equilibrium lattice constant ao (Tab.QJ. Site- 
decomposed densities for X, B, Sc are plotted in brown, red 
and blue, respectively. The Fermi level is shifted to zero and 
marked by a vertical line. 



E(k) along the X-M-T direction. Generally, the bands 
in ScglnB are more dispersive comparing e.g. to MgCNia 
(53|,j5^ or Scain (see Fig.^ll. Two separated, lowest ly- 
ing bands are formed from s-states of X and B atoms, 
with notably large contribution from Sc. 

For our searching, the most interesting feature of 
electronic structure is related to quite large McMillan- 



Atom 


n{EF) 








A, 


Sc 


20.0 


18 


15.2 


50 


0.36 


Tl 


17.4 


2 


15.2 


226 


0.01 


B 


7.6 


32 


61.9 


204 


0.16 


Sc 


20.0 


19 


18.5 


76 


0.25 


In 


16.3 


1 


18.5 


193 


0.01 


B 


6.7 


34 


58.3 


181 


0.19 


Sc 


19.0 


18 


19.0 


80 


0.23 


Ga 


18.4 


2 


19.0 


124 


0.01 


B 


6.5 


35 


60.8 


197 


0.18 


Sc 


21.6 


18 


21.3 


100 


0.18 


Al 


21.5 


1 


21.3 


60 


0.02 


B 


6.4 


30 


58.6 


183 


0.16 



Hopfield parameters seen along the whole series 
(see Tab. HH) . In scandium, typically for transition ele- 
ments, d - f scattering channel gives the most important 
contribution to rji. For boron and X-element only the 
p - d channel contributes to Note that although the 
rjB values were found to be the largest in our systems, 
•qsc occurred to be more important in calculation of A, 
since it is counted three times (3 Sc atoms in unit cell), 
and has lower value of Mi{ujf) (see below), rjx has neg- 
ligible value, despite a noticeable density of states at the 
Fermi level. The accuracy of computed rji, associated 
with different MT geometry, is about 1 mRy/Bohr^ for 
scandium and 2 mRy/Bohr^ for boron. 

The phonon DOS, and its evolution with the X- 
element, is given in Fig. 01 The most striking feature 
of the presented F{uj) is the 'rigid-band-like' modifica- 
tion with the mass variation of the X-element (note that 
the mass distribution in the unit cell markedly changes) . 
The highest peak in Fig. 0] marks the lowly-dispersive 
part of acoustic branches, associated with the X-atom 
vibrations. As one can see, when the mass of X is get- 
ting smaller, the range of acoustic modes broadens, which 
is manifested in a shift of their flat parts towards higher 
frequencies, i.e. the peak moves from 9 meV {X=T\) 
to 27 meV {X=A\). Generally, the phonon spectrum 
consists of two separate areas: the high-frequency part 
includes essentially phonon modes of light B, with en- 
ergies above 50 meV, and the lower-frequency part, be- 
low 50 meV, contains mainly X and Sc states. This be- 
havior was also characteristic for the dynamic properties 
of MgCNia compound, where carbon vibrations had the 
highest frequencies ^55,]. 

Let us remark, that in the phonon calculations for all 
compounds we met some problems with imaginary fre- 
quencies, appearing near the BZ center (F), and addition- 
ally for X — T\ and Ga, near R point (in Fig. 0] they are 
visible as negative frequencies tails). These eigenvalues 
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FIG. 3: Electronic dispersion curves E(k) along high symme- 
try directions in ScalnB perovskite. The Fermi level is shifted 
to zero and marked by a horizontal line. 



correspond to three optical modes, and at F and R points 
only Sc atoms vibrate in these modes, with eigenvectors 
lying in three planes perpendicular to the Sc-B bonds. 
Physically, the occurrence of imaginary frequencies may 
indicate either the instability of the perfect perovskite 
structure (favoring crystal distortion) or important an- 
harmonic effects. 

Similar^, unstable phonon branches were detected in 
MgCNia jSa, [s^l occurring at different high-symmetry 
points (X, M). These branches corresponded mainly to 
anharmonic Ni vibrations, resulted from a double-well 
potential [s^ [s^ [s^, in which Ni was placed. How- 
ever, due to the shallowness of the double well, no stable 
long-range structural distortion was found there, and the 
perovskite structure was stabilized dynamically. The ori- 
gin of unstable vibrations in SC3XB is attractive problem 
itself, but more detailed analysis of their dynamical prop- 
erties is out of the scope of this paper. 

In summary, the influence of the negative frequency 
range on phonon density of states is not so critical (the re- 
lated energy bands are strongly dispersive). This should 
not affect much the (cuf) calculations. The total weight 
of the negative frequency area is about 0.5% in ScalnB 
and SC3AIB, 2.5% in ScgGaB and 5% in SC3TIB. So, in 
the case of X = In and Al the 'negative tail' was negli- 
gible and cut off in the frequency moments calculations. 
For two remaining compounds {X — Ga and Tl) the tail 
was cut off at about 3 meV and F{uj) was extrapolated to 
reach a; = in Debye-like manner. Because the phonon 
DOS obtained for SC3TIB seems less reliable, in this case 
the evaluation of the superconducting parameters should 
be treated rather qualitatively. 

Having phonon DOS we may proceed to estimate EPC 
strength. The mean phonon frequencies for constituent 
atoms were computed from the phonon DOS using the 
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FIG. 4: Phonon DOS F{uj) in ScaXB perovskites. 



above-mentioned analysis of the spectrum: (uif) for X 
and Sc were calculated from the lower-frequency part of 
the spectrum, and were assumed to be the same for both 
atoms. As expected, the lowest value of (ujf) is observed 
in thallium compound. The average (ujg) (for boron) 
were computed from the high-frequency part of DOS, 
and resulted in significantly larger values, due to small 
mass of B (see Tab. Hlj) . Note, that taking the same (ujf) 
for all atoms (sometimes practiced), might occur incor- 
rect [31I I and disagree with atomic-like character of r]i 
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TABLE III: Total EPC constant A and critical temperature 
(K) for fi* = 0.13. (lu) used in Eq. @ is given in (K). 





SC3TIB 


ScglnB 


ScgGaB 


SC3AIB 


A 


1.25 


0.94 


0.88 


0.72 


(lo) 


182 


225 


226 


263 




15 


12 


10 


7.5 



TABLE IV: Influence of unit cell volume shrinking on elec- 
tronic critical parameters. 



Compound 


dn{Ep) 






dV 


rfln V 


dh^V 


SC3TIB 


0.41 


-1.4 


-1.9 


Sc3lnB 


0.41 


-1.5 


-2.0 


Sc3GaB 


0.50 


-1.2 


-2.3 


SC3AIB 


0.44 


-1.4 


-2.6 



{MiLof corresponds to an effective force constant, so it 
is also a site-dependent quantity). In systems, where 
atoms fiave markedly different masses, tfie lightest el- 
ement would be favored in that way, as shown below. 
Therefore, the phonon DOS was divided into high- and 
low- frequency parts in the {ujD calculations. 

The values of total EPC constant A, i.e. sum over 
atomic contributions shown in Tab. ^ are gathered in 
Tab. mil together with the critical temperature values. 
In all investigated compounds the main contribution to 
A comes from scandium. Due to the lowest value of (1^5^,) 
parameter, SC3TIB has the highest A. Superconducting 
transition temperature Tc was estimated using Eq. 
with typical value of Coulomb pseudopotential ^* — 0.13 
applied. As one can notice, all compounds with calcu- 
lated A in the range of 0.7 - 1.2, and Tc ~ 10 K are 
medium- or even strong-coupling superconductors in the 
RMTA framework. 

To verify the influence of employed //* on Tc magni- 
tude, we plotted the value of Tc versus /i* in the rea- 
sonable range of 0.09 < < 0.21 (Fig. Transi- 
tion temperatures decrease almost linearly with fj,* in 
this range, and even for large value of fi* = 0.21, Tc 
is still enough high to be detectable in the typical low- 
temperature measurements. 
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Sc,AlB 
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FIG. 5: Influence of /i* parameter on critical temperature 
Tc. Vertical line marks fi* = 0.13. 

We also examined how our results changed, when (utf) 
was taken to be equal for all atoms. For example, in 
ScalnB we obtaine d A — 1.6 with As = 1.1 and Xsc — 0.5, 
when employed ^/ (uj^) — 23.4 meV, as derived from the 
entire phonon spectrum. The resulting critical tempera- 
ture Tc = 26 K was much higher than the previous value, 
and probably much overestimated, due to the overesti- 



mation of boron contribution. 

The pressure effect on electronic properties of SC3XB 
was additionally inspected by decreasing lattice param- 
eter in KKR computations (in the range of 0-5%). The 
variation of the most important parameters - total den- 
sity of states at Ep and McMillan-Hopfield parameters 
for boron and scandium - is summarized in Tab. lIVl The 
observed tendency is again similar in the whole series. In 
spite of the n[Ep) decrease with volume shrinking, rji pa- 
rameters increase with ratio about dlnrjB/dlnV ~ —2 
and dlnrjsc/ d\nV ~ —1.5. 

These results suggest, that external pressure may en- 
hance Tc, as observed in MgCNia 0,0, but the opposite 
effect of lattice stiffening may be dominant. 

To have a better insight into the magnitude of 
McMillan-Hopfield parameters in SC3XB, we performed 
similar electronic structure calculations for the exist- 
ing, related perovskite superconductors - MgCNia and 
YRh3B, applying the experimental values of lattice pa- 
rameters (in Bohr units): 7.205 (Ref. 0) and 7.870 
(Ref. 1^), respectively. Since DOS for MgCNis was re- 
cently published by many authors (e.g. Refs. 53, 5^1), we 
present KKR DOS for the less known YRhsB in Fig. El 



TABLE V: McMillan-Hopfield parameters for MgCNi3 and 
YRh3B (mRy/BohrVatom). 



Element 


V 


Element 


V 


Ni 


20.4 


Rh 


17.5 


Mg 


0.1 


Y 


1.4 


C 


9.1 


B 


1.6 



As we can observe from Tab. ^ and Tab. values 
of rji for the transition element (Sc, Ni, Rh) are very 
similar in all cases. It is interesting to underline that 
the presence of the 2p-states of light element (B in this 
case) at the Fermi level in YRh3B is negligible [5ll IH^ ] 
(see Fig. EJ, unlike in MgCNi3 or SC3YB. This is the 
reason of a very small rjB parameter, comparing to rjc 
in MgCNi3 and especially rjB in SC3YB. Also the large 
mass of rhodium [M ~ 103 u), comparing to scandium 
{M ~ 45 u) or nickel [M ~ 59 u) may be responsible for 
the low superconducting critical temperature (~1 K) in 
this compound. The comparison of McMillan-Hopfield 
parameters in presented perovskites additionally favors 
superconductivity in ScsATB compounds due to relatively 
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FIG. 6: KKR DOS of YRhaB calculated at experimental 
lattice parameter. Site-decomposed densities for Y, B, Rh 
are plotted in brown, red and blue, respectively. The Fermi 
level is shifted to zero and marked by a vertical line. 



low mass of scandium. 



B. Effect of vacancy on B-site 

When dealing with this family of compounds, one has 
to face with possible crystallographic imperfections of the 
SC3XB structure. The common problem occurring in in- 
termetallic perovskite borides (and carbides) is related 
to not fully occupied boron (carbon) position. Note that 
various boron or carbon deficiency was observed in many 
related perovskites '^^ll as in MgCNis. This 

effect is very important for superconductivity, since in 
MgC^Nis critical temperature decreases linearly with C 
concentration 58] and superconductivity disappeares for 
X < 0.9. The vacancy on B-site is highly plausible in our 
compounds and this effect has to be taken into account 
in the analysis. This prompted us to perform the cal- 
culations of electronic structure and McMillan-Hopfield 
parameters in ScaXB^: from the KKR-CPA method [s^ . 
A vacancy on boron site was treated as an 'empty sphere' 
with Z = 0, and the same MT-radius as applied for B 
atom. In the non-stoichiometric ScaXB^; the (lb) posi- 
tion (1,5,5) is occupied by B atom and a vacancy with 
probabilities x and 1 - a;, respectively. The electronic 
structure of such disordered system was calculated using 
the coherent potential approximation (CPA), which al- 
lows considering any finite concentrations, including im- 
purity states. 

Site-decomposed DOS at Ep in Sc^XB^ (divided by 
nl^Ep) in stoichiometric SC3XB) and McMillan-Hopfield 
parameters, versus B concentration, are shown in Fig. 
At a first glance one observes, that all site contributions 
to DOS at Ep decrease in the investigated systems with 



decreasing of boron concentration, and reach less than 
50% of initial values for x = 0.85. This is unlikely to the 
rigid-band behavior, where Ep is expected to move to 
the left, towards higher DOS (see Fig. O. More detailed 
analysis of KKR-CPA DOS in ScgXB^ (see also Fig. [El) 
indicates, that a vacancy on B-site seems to behave as a 
hole donor, if the pseudo-gap in DOS, found above Ep, 
can be considered as a separation between valence-like 
and conduction-like bands. Since the potential of the 
vacancy is much more repulsive than the potential of B 
atom, all p-states (accommodating six electrons) are ex- 
pelled into higher energy range (well above Ep), against 
only one electron occupying p-shell in B atom. Con- 
sequently, the filling of low-lying conduction-like states 
decreases, when the vacancy concentration increases. 

The critical parameters analysis for considered vacancy 
concentrations (0.85 < a; < 1) shows that presence of 
vacancy is very unfavorable for occurrence of supercon- 
ductivity in these structures, due to a sudden decrease 
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spectively. Lines connecting calculated points are added as a 
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of the most important rjsc and t]b parameters (Fig. [T]). 
Using the values of {ojf) obtained in the stoichiomet- 
ric ScsXB compounds, the variation of total A and Tc 
for boron-deficient structures was analyzed. As we can 
see in Fig. |S| in all structures the decrease in the EPC 
constant A is so fast, that when vacancy concentration 
reaches 15% (i.e. SC3XB0.85) superconductivity is prac- 
tically quenched, since Tc ^ 0.1 K. So, the boron atom 
occurs to play crucial role in superconductivity onset in 
SC3XB, and even small boron deficiency may cause rapid 
decrease of the critical temperature. 

This behavior seems to be similar to MgCNia case, 
where lowering of Tc with increase of vacancy concentra- 
tion on carbon site, was observed experimentally 58, 5^. 
Quite recently, two important investigations were car- 
ried for this system. It was deduced from the specific 
heat measurements |60| that in a carbon-deficient sam- 
ple low-energy nickel phonon modes (which are probably 
the most important in superconductivity) were shifted to 
higher energies. Furthermore, electronic structure calcu- 
lations performed for MgC^^Nia in the range 0.80 < x < 1 
exhibited the decrease of n{Ep) with the increasing x 
[6ll |. but less drastic than in our ScaXB^: systems. So, 
disappearance of superconductivity in MgCajNia seems to 
result from two unfavorable effects: lowering the DOS at 
Fermi level combined with increasing phonon frequencies. 
Taking into account the possibility of similar behavior of 
low-energy phonon modes in Sc^X'Qx, we may suggest, 
that superconductivity in SC3XB could be even more sen- 
sitive to stoichiometry, than in MgCNis. 



Compound 


a experimental [16] 


ao calculated 


SC3TI 




8.300 


Scsin 


8.427 


8.150 


ScsGa 




8.025 


SC3AI 




8.150 



C. Magnetic properties of cubic SC3X 

The motivation to study relations between supercon- 
ductivity in SC3XB, and magnetism in SC3X compounds, 
was inspired by the widely-studied weak itinerant ferro- 
magnetism of the hexagonal Scain . In spite of recent 
interest in this field, the research of the cubic form of 
ScaIn has not been carried out so far. To have a possibil- 
ity of wider comparison between the two series of com- 
pounds: ScaXB and ScaX, the electronic structure cal- 
culations for three hypothetical structures (SC3AI, ScaGa 
and ScaTl) were also performed. 

In all cubic Sc^X structures lattice constants were de- 
rived from the total energy calculations, and MT sphere 
radii for both Sc and X atoms were set to 0.35ao (other 
computational details are the same as in Sec. Ill All . Ta- 
ble El presents the calculated values of equilibrium lat- 
tice parameter (ao). The variation of oq with X-element 
in SC3X is similar to the tendency observed in Sc3ArB. 
However, there is much larger difference (3%) between 
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FIG. 9: KKR non-polarized total DOS in CusAu-type 
ScaGa and SC3AI calculated at equilibrium lattice constant 
(Tab.Elll. The Fermi level is shifted to zero and marked by 
a vertical line. 
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FIG. 10: KKR polarized total DOS in CusAu-type ScsTl and 
Scsin calculated at equilibrium lattice constant (Tab. IVI^ . 
The Fermi level is shifted to zero and marked by a vertical 
line. 



the theoretical and experimental values for Scain, but 
still within the acceptable error of LDA. This however 
prompted us to study electronic structure as a function 
of lattice constant in this compound. 

First, we have analyzed the non-polarized DOS and 
the Stoner product for SC3X compounds (Tab. IVII|I . 
In view of the KKR results, the Stoner criterion is 
not fulfilled only in ScaGa and SC3AI, and indeed, the 
spin-polarized KKR computations yielded non-magnetic 
ground state in both cases (non-spin-polarized DOS are 
shown in Fig. O. 

In agreement with the Stoner analysis, spin-polarized 
KKR computations in SC3TI and Scain converged to fer- 
romagnetic ground state (Fig. I10|l with magnetic mo- 
ments as large as 0.37 /is and 0.25 /is per scandium atom 
for the case oi X — T\ and In, respectively. The magnetic 
moment on X atom was negligible (~ 0.02 /is). Regard- 
ing to the calculations, scandium in cubic phase of Scain 
possesses quite similar magnetic moment as in the hexag- 
onal one [73: HSc = 0.32 /i^ (comparing with other LDA 



TABLE VIL Electronic properties of SciX. n{EF) is 
given in (1/Ry/spin/f.u.) (f.u. - formula unit), nsc in 
(1/Ry/spin/atom), nsc in (/.is/atom). 



non-spin-polarized calculations 


Compound 


n{EF 


) 


nsc{EF) 


/ X n{EF) 


SC3TI 


91.2 




30.1 


1.78 


Scgln 


62.1 




19.9 


1.17 


ScsGa 


45.9 




14.4 


0.89 


SC3AI 


37.2 




11.8 


0.70 


Scsin (hex) 


72.8 




23.0 


1.32 


spin-polarized calculations 


Compound 


n|(_EF) n 


(Ef) 


wtSc(-Bf) nisciEp) 




SC3TI 


43.3 


26.4 


14.2 8.5 


0.37 


Scgln 


49.6 


34.7 


15.7 11.0 


0.25 


Scsin (hex) 


69.7 


27.5 


19.3 8.5 


0.32 



result 0.30 hb ^m)- The magnitude of magnetic mo- 
ment in Scsin as a function of lattice constant is shown 
in Fig 111! The observed monotonic decrease of magne- 
tization with the cell volume shrinking indicates, that 
the quantum critical point, when magnetization disap- 
pears, can be reached under hydrostatic pressure. Con- 
versely, magnetization in hexagonal ScaIn enhances [b^ 
under moderate hydrostatic pressurCjWhich is supported 
by electronic structure calculations [63 suggesting an in- 
crease of magnetic moment under weak hydrostatic pres- 
sure, while suppression of ferromagnetism upon applying 
the uniaxial strain. 

It seems interesting to recall, that the widely inves- 
tigated weak ferromagnet NiaAl ■64] crystallizes in the 
same CuaAu-type structure. Also, the electronic struc- 
ture calculations in NiaAl gave weak magnetic moments 
on Ni {fiNi — 0.24 fj,B 113), which is close to the KKR 
value we have gained for scandium in the cubic ScaIn 
(/iSc — 0.25 hb)- According to measurements [i^, NiaAl 
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FIG. 11: Magnetic moment on scandium atom in Scsin as a 
function of the lattice constant. Line is added as a guide for 
the eye. 
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exhibits weak itinerant magnetism, with however much 
smaller magnetic moment fi^i = 0.075 ^b- Surprisingly, 
the similar compound NiaGa was found to be paramag- 
netic in experiment [g^I, although LDA calculations re- 
sulted in the magnetic state with /i^vi = 0.26 /is 65\, even 
larger than in NisAl. The discrepancy between LDA pre- 
diction and experimental finding in NisGa was explained 
in terms of strong spin fluctuations 65] . 

The question whether ScsX systems are similar to 
NiaA" compounds will be interesting after experimental 
verification of existence of these systems. 



D. Ferromagnetism vs. Superconductivity 

As we can see from the limiting cases, studying possi- 
ble competition of superconductivity and weak ferromag- 
netism is directly related to investigation of the role of 
trivalent boron in the entitled compounds. To enlighten 
this subject, the KKR-CPA calculations were performed 
in full range of B concentration for illustrative example 
of ScglnB 

Analyzing the non-magnetic DOS values for Scain 
(Tab. IVIII) one notices, that the large n{Ep) value con- 
sists mainly (over 95%) of scandium atoms contribution. 
This gives large Stoner product on Sc and makes the 
magnetic ground state energetically favorable. The DOS 
evolution (Fig. I13|l shows that when B concentration in- 
creases, the Fermi level shifts from the strongly increas- 
ing DOS (5% B in Fig.^J, leaving the 'magnetic' region, 
towards a deep valley (75% B in Fig. [T^ . After cross- 
ing this minimum, the n{Ep) value again increases, but 
the Stoner limit on Sc is not reached here (95% B in 
Fig. I13|l . The total DOS at Ep, as well as the scandium 
contribution nsc{Ep) in ScalnB is much lower, than the 
corresponding values in Scsln. Consequently, the Stoner 
product in ScalnB is as small as 0.5, which prevents for- 
mation of magnetic ground state. 

In order to better understand mechanism of lowering 
the scandium contribution to n{Ep), we can notice inter- 
esting result concerning the formation of one additional 
low-lying band in ScalnB (two lowest-lying bands below 
-0.3 Ry in Fig. |3| against only one in Scain in Fig. [T^ . 
This band is formed from hybridization of s— like states 
from B with Sc orbitals (including d— states) and presum- 
ably bounds about one electron from B and one electron 
from Sc. Consequently, the energy bands in ScalnB in 
the range of —0.3 Ry < E < Ep) includes one Sc elec- 
tron less, than the bands in Scain in the corresponding 
energy range. 

When comparing integrated partial c?-DOS for Sc atom 
in both compounds, we clearly notice the decrease of 
d— Sc orbitals filling in ScalnB: 1.5 e against 1.7 e in 
Scain (per Sc atom) as well as the transfer of about 0.1 e 
from upper bands to the lower additional band for each 
Sc atom. 

These two effects show the important modification of 
valence states of Sc upon B insertion and seem also to 
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FIG. 12; The non-polarized electronic dispersion curves E(k) 
along high symmetry directions in cubic Scsin (for ao = 8.150 
Bohr). The Fermi level is marked by a horizontal line. 



be responsible for lowering the scandium contribution to 
n{Ep) and preventing the ferromagnetism. When the 
ferromagnetic ground state in ScaTlB and ScalnB is de- 
stroyed, the conventional superconductivity is allowed to 
appear. Certainly, the presence of boron atom is not 
only important for preventing ferromagnetism, but, as 
we have seen in Sec. Ill Bl it is crucial to promote super- 
conductivity. This opens the perspectives of interesting 
experimental study, if it would be possible to synthesize 
the ScaXBa: systems with various boron concentration. 

Appearance of magnetic ground state in ScaTl and 
ScaIn makes ScaTlB and ScalnB even more similar 
to MgCNia, since hypothetical CuaAu-type structure 
MgNia was theoretically predicted to have a small mag- 
netic moment (about 0.4 /is per Ni atom l6^). We may 
confirm this result {fj,Ni = 0.40 fj,B from KKR). 

For comparison with ScalnB^; we present the evolu- 
tion of electronic structure for MgC^jNia fFig. I13|l . The 
DOS variation with C concentration is a bit different, 
since Ep moves towards lower lying valence states when 
the carbon concentration decreases, in contrast to the 
ScalnBa; system. However, the general trends are quite 
similar in both cases. When C concentration decreases 
the Fermi level crosses the DOS valley and next falls into 
the higher-DOS region, where magnetic ground state ap- 
pears. Also it is worth noting, that van Hove singularity 
near Ep in MgCNia, being the unusual feature of its elec- 
tronic structure, disappears with the vacancy concentra- 
tion (due to the deficiency of carbon p-orbitals). Note 
that our DOS picture for x — 75% compares well with 
the result obtained for x — 80% in Ref. |61| . 

Interestingly, the comparison of data presented in 
Tab, mil and Tab. IVIll suggests. that the strength of mag- 
netic interactions in ScaA^ and the electron-phonon in- 
teractions in ScaATB seems to be correlated, since both 
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FIG. 13: Effect of vacancy on total density of states in ScsInB^; (upper panels) and MgC^Nis (lower panels) from KKR-CPA 
calculations. Note a remarkable evolution of electronic states in the vicinity of the Fermi level: from high DOS at Ef in 
C-rich and B-rich samples (left side: i.e. near full C/B occupancy and superconductivity limit) through a deep DOS minimum 
(intermediate concentrations, i.e. near 50% C in MgCaiNis and c.a. 80% B in ScsInB^;) to again large DOS, which satisfy the 
Stoner criterion (right side, i.e. near empty C/B sublattice limit). The Fermi level is at zero, marked by a vertical line. 



the Stoner product / x n^Ep) and EPC parameter A are 
increasing with X, i.e. the lowest values are observed for 
X = A\ and highest for X = Tl. 



III. CONCLUSIONS 

We have presented theoretical investigation of super- 
conducting properties of the perovskite series SC3XB 
{X = Tl, In, Ga, Al) and their possible connections with 
weak magnetism in the corresponding SC3X compounds. 
Our main results predicted: 

(i) superconductivity in SC3XB, with A ~ 0.7 - 1.25 and 
Tc ~ 7 - 15 K; 

(ii) weak ferromagnetism in Scain and SC3TI and the ab- 
sence of ferromagnetism in Sc3Ga and SC3AI; 

(iii) critical effect of vacancy on B-site on superconduc- 
tivity in Sc3XBa;. 

We have also showed that boron inserted to the cubic 
SC3TI and Scsin destroys the magnetic ground state, and 
likely turns these systems into superconductors. 

On the whole, the estimated EPC parameters and criti- 
cal temperatures Tc obtained for SC3XB compounds from 
the RMTA analysis are even larger that the values re- 
ported for two perovskite superconductors MgCNi3 and 
YRh3B. 



Preliminary experimental examination of one of the 
compounds from the entitled series - ScsInB - was un- 
dertaken but several synthesis procedures did not 
succeed in preparing a single-phase and stoichiometric 
compound. However, one of the samples showed signs of 
superconductivity with Tc ~ 4.5 K. Further experimental 
study is needed to clarify this result. 

Moreover, in view of the fact that some itinerant ferro- 
magnets exhibit superconducting properties, e.g. Y9C07 
[67| or UGe2 [i^ , the experimental investigation of cubic 
series SC3X should also be very appealing. 



12 



[1] T. He, Q. Huang, A. P. Ramirez, Y. Wang, K. A. Regan, 

N. Rogado, M. A. Hayward, M. K. Haas, J. S. Slusky, 

K. Inumara, et al.. Nature 411, 54 (2001). 
[2] H. Rosner, R. Weht, M. D. Johannes, W. E. Pickett, and 

E. Tosatti, Phys. Rev. Lett. 88, 27001 (2002). 
[3] T. G. Kumary, J. Janaki, A. Mani, S. M. Jaya, V. S. 

Sastry, Y. Hariharan, T. S. Radhakrishnan, and M. C. 

Valsakumar, Phys. Rev. B 66, 64510 (2002). 
[4] G. Garbarino, M. Monteverde, M. Nez-Regueiro, 

C. Acha, R. Weht, T. He, K. A. Regan, N. Rogado, 

M. Hayward, and R. J. Cava, Physica C 408, 754 (2004). 
[5] Z. Q. Mao, M. M. Rosario, K. D. Nelson, K. Wu, I. G. 

Deac, P. SchiflFer, Y. Liu, T. He, K. A. Regan, and R. J. 

Cava, Phys. Rev. B 67, 94502 (2003). 
[6] R. Prozorov, A. Snezhko, T. He, and R. J. Cava, Phys. 

Rev. B 68, 180502 (2003). 
[7] J. Y. Lin, P. L. Ho, H. L. Huang, P. H. Lin, Y. L. Zhang, 

R. C. Yu, C. Q. Jin, and H. D. Yang, Phys. Rev. B 67, 

52501 (2003). 

[8] P. M. Singer, T. Imai, T. He, M. A. Hayward, and R. J. 

Cava, Phys. Rev. Lett. 87, 257601 (2001). 
[9] T. Klimczuk and R. J. Cava, Phys. Rev. B 70, 212514(R) 

(2004). 

[10] H. Takei, N. Kobayashi, H. Yamauchi, T. Shishido, and 
T. Fukase, J. Less-Common Met. 125, 233 (1986). 

[11] R. E. Schaak, M. Avdeev, W.-L. Lee, G. Lawes, H. W. 
Zandbergen, J. D. Jorgensen, N. P. Ong, A. P. Ramirez, 
and R. J. Cava, J. Solid State Chem. 177, 1244 (2004). 

[12] M.-S. Park, J. Giim, S.-H. Park, Y. W. Lee, S. L Lee, 
and E. J. Choi, Supercond. Sci. Technol. 17, 274 (2004). 

[13] M. D. Johannes and W. E. Pickett, Phys. Rev. B 70, 
060507(R) (2004). 

[14] B. Wiendlocha, J. Tobola, S. Kaprzyk, D. Fruchart, 
and J. Marcus, The European Conference Physics of 
Magnetism'05, 24-27.06.2005, Poznan, Poland. Preprint: 
cond-mat/0506793 

[15] H. HoUeck, J. Less-Common Met. 52, 167 (1977). 

[16] A. Palenzona, P. Manfrinetti, and R. Palenzona, J. All. 
Compd. 15, 94 (1996). 

[17] J. F. Cannon, High Pressure in Science and Technol- 
ogy. Proceedings of the 9-th AIRAPT International High 
Pressure Conference p. 113 (1984). 

[18] B. T. Matthias, A. M. Clogston, H. J. Wilhams, 
E. Corenzwit, and R. C. Sherwood, Phys. Rev. Lett. 7, 
7 (1961). 

[19] W. L. McMillan, Phys. Rev. 167, 331 (1968). 
[20] J. J. Hopfield, Phys. Rev. 186, 443 (1969). 
[21] G. D. Gaspari and B. L. Gyorffy, Phys. Rev. Lett. 28, 
801 (1972). 

[22] W. E. Pickett, Phys. Rev. B 25, 745 (1982). 

[23] L L Mazin, S. N. Rashkeev, and S. Y. Savrasov, Phys. 
Rev. B 42, 366 (1990). 

[24] A. D. Zdetsis, E. N. Economou, and D. A. Papaconstan- 
topoulos, Phys. Rev. B 24, 3115 (1981). 

[25] W. H. Butler, J. J. Olson, J. S. Faulkner, and B. L. Gy- 
orffy Phys. Rev. B 14, 3823 (1976). 

[26] W. H. Butler, Phys. Rev. B 15, 5267 (1977). 

[27] S. Y. Savrasov and D. Y. Savrasov, Phys. Rev. B 54, 
16487 (1996). 

[28] D. A. Papaconstantopoulos, L. L. Boyer, B. M. Klein, 
A. R. Williams, V. L. Morruzzi, and J. F. Janak, Phys. 



Rev. B. 15, 4221 (1977). 

[29] S. S. Rajput, R. Prasad, R. M. Singru, S.Kaprzyk, and 
A. Bansil, J. Phys.: Condens. Matter 8, 2929 (1996). 

[30] B. M. Klein, L. L. Boyer, and D. A. Papaconstantopou- 
los, Phys. Rev. Lett. 42, 530 (1979). 

[31] B. M. Klein and D. A. Papaconstantopoulos, Phys. Rev. 
Lett. 32, 1193 (1974). 

[32] S. V. Halilov, D. J. Singh, and D. A. Papaconstantopou- 
los, Phys. Rev. B 65, 174519 (2002). 

[33] W. E. Pickett and D. J. Singh, Phys. Rev. Lett. 72, 3702 
(1994). 

[34] D. A. Papaconstantopoulos and B. M. Klein, Phys. Rev. 

Lett. 35, 110 (1975). 
[35] J. Kortus, L I. Mazin, K. D. Belashchenko, V. P. 

Antropov, and L. L. Boyer, Phys. Rev. Lett. 86, 4656 

(2001). 

[36] W. E. Pickett, Rev. Mod. Phys. 61, 433 (1989). 

[37] H. J. Choi, D. Roundy, H. Sun, M. L. Cohen, and S. G. 

Louie, Nature 418, 758 (2002). 
[38] A. Floris, G. Profeta, N. N. Lathiotakis, M. Luders, 

M. A. L. Marques, C. Franchini, E. K. U. Gross, A. Con- 

tinenza, and S. Massidda, Phys. Rev. Lett. 94, 37004 

(2005). 

[39] A. Bansil, S. Kaprzyk, P. E. Mijnarends, and J. Tobola, 

Phys. Rev. B 60, 13396 (1999). 
[40] W. H. Butler, P. Dederichs, A. Gonis, and R. Weaver, 

eds.. Chapter HI, in: Applications of Multiple Scattering 

Theory to Materials Science, vol. 253 (MRS Symposia 

Proceedings, MRS Pittsburgh, 1992). 
[41] S. B. der Kellen, Y. Oh, E. Badralexe, and A. J. Freeman, 

Phys. Rev. B 51, 9560 (1995). 
[42] M. Ogura and H. Akai, J. Phys.: Condens. Matter 17, 

5741 (2005). 

[43] T. Stopa, S. Kaprzyk, and J. Tobola, J. Phys.: Condens. 

Matter 16, 4921 (2004). 
[44] J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 

(1992). 

[45] S. Baroni, A. D. Corso, S. de Gironcoli, and P. Giannozzi, 

www . pwscf . org . 
[46] S. Baroni, A. D. Corso, S. de Gironcoli, and P. Giannozzi, 

Rev. Mod. Phys 73, 515 (2001). 
[47] J. P. Perdew and A. Zunger, Phys. Rev. B 23, 5048 

(1981). 

[48] M. Methfessel and A. T. Paxton, Phys. Rev. B 40, 3616 
(1989). 

[49] P. B. Allen and R. C. Dynes, Phys. Rev. B 12, 905 (1975). 
[50] G. Grimvall, The electron-phonon interaction in metals 

(North-Holland Publishing CO. Amsterdam, 1981). 
[51] P. Ravindran, S. Sankaralingam, and R. Asokamani, 

Phys. Rev. B 52, 12921 (1995). 
[52] P. J. T. Joseph and P. P. Singh, cond-mat/0505310 vl. 
[53] S. B. Dugdale and T. Jarlborg, Phys. Rev. B 64, 

100508(R) (2001). 
[54] A. Szajek, J. Phys.: Condens. Matter 13, 595 (2001). 
[55] R. Held, B. Renker, H. Schober, P. Adelmann, D. Ernst, 

and K. P. Bohnen, Phys. Rev. B 69, 92511 (2004). 
[56] A. Y. Ignatov, S. Y. Savrasov, and T. A. Tyson, Phys. 

Rev. B 68, 220504 (2003). 
[57] A. Y. Ignatov, L. Dieng, T. Tyson, T.He, and R. Cava, 

Phys. Rev. B 67, 064509 (2003). 
[58] T. Amos, Q. Huang, J. Lynn, T. He, and R. Cava, Solid 



13 



State Comm. 121, 73 (2002). 
[59] Z. A. Ren, G. C. Che, S. L. Jia, H. Chen, Y. M. Ni, G. D. 

Liu, and Z. X. Zhao, Physica C 371, 1 (2002). 
[60] A. Wahe, G. Fuchs, K.-H. MuUer, S.-L. Drechsler, 

K. Nenkov, and L. Schuhz, Phys. Rev. B 72, 100503(R) 

(2005). 

[61] P. J. T. Joseph and P. P. Singh, Phys. Rev. B 72, 064519 
(2005). 

[62] A. Aguayo and D. J. Singh, Phys. Rev. B 66, 20401(R) 
(2002). 

[63] J. Grewe, J. S. Schilhng, K. Ikeda, and J. K. A. Gschnei- 

dner, Phys. Rev. B 40, 9017 (1989). 
[64] F. R. de Boer, C. J. Schinkel, J. Biesterbos, and 

S. Proost, J. Appl. Phys. 40, 1049 (1969). 
[65] A. Aguayo, I. I. Mazin, and D. J. Singh, Phys. Rev. Lett. 

92, 147201 (2004). 
[66] J. H. Shim, S. K. Kwon, and B. I. Min, Phys. Rev. B 64, 

180510 (2001). 

[67] A. Kolodziejczyk and J. Spaiek, J. Phys. F 14, 1277 
(1984). 

[68] S. S. Saxena, P. Agarwal, K. Ahilan, F. M. Grosche, 
R. K. W. Hasselwimmer, M. J. Steiner, E. Pugh, I. R. 
Walker, S. R. Juhan, P. Monthoux, et al.. Nature 406, 
587 (2000). 

[69] V. B. Compton and B. T. Matthias, Acta Cryst. 15, 94 
(1962). 

[70] For a review of the MgCNis properties see also the recent 



papers on this subject, e.g. Refs. |gll6l|. 

[71] Definition of the n-th frequency moment of the Eliash- 
berg coupling function a'^F{Lo) [i^, [H^j: (t^^") = 
/ Lu''-^a'^F{uj)duj/ J Lu~^a^F{Lo)dLu. If a^{uj) ~ const, 
{lo") ~ / Lo"~^F{ij)duj/ J Lj~-^F{u)dLj. a^{u) is the 
electron-phonon interaction coefficient. 

[72] The values of n(EF) presented in Tab.Ellfor ScsInB are a 
bit different from those in Ref . [T^ . It comes from applied 
smaller lattice parameter (the equilibrium value instead 
of the experimental one) and semi-relativistic calcula- 
tions. Because the Fermi level is located on the slope of 
DOS peak, small shift in Ef has noticeable influence on 
n{EF). However, rj values are not infected much, because 
they are defined as a ratio of resulting densities. 

[73] In the NisSn-type Scsin experimental lattice constants 
were taken: a — 6.42 A, c = 5.18 A |69|l . the adjustable 
parameter x, positioning Sc atoms in the XY plane, was 
assumed to have the ideal value of j: = 5/6, following the 
X-ray analysis |69| and previous calculations i62l . The 
same MT sphere radius (3.0 Bohr) were used for Sc and 
In. In the present work the magnetic moment is given per 
Wigner-Seitz cell and per Sc atom. In Ref. [T^ the mag- 
netic moments of 0.26 hb (hexagonal phase) and 0.27 /is 
(cubic phase, experimental lattice constant), were given 
per Sc muffin-tin sphere. 



